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Abstract. We study the Schrodinger equation of a class of two-level systems under the action 
of a periodic time- dependent external field in the situation where the energy difference 2e between 
the free energy levels is sufficiently small with respect to the strength of the external interaction. 
Under suitable conditions we show that this equation has a solution in terms of converging power 
series expansions in e. In contrast to other expansion methods, like in the Dyson expansion, the 
method we present is not plagued by the presence of "secular terms". Due to this feature we were 
able to prove absolute and uniform convergence of the Fourier series involved in the computation 
of the wave functions and to prove absolute convergence of the e-expansions leading to the "secular 
frequency" and to the coefficients of the Fourier expansion of the wave function. 
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1 Introduction 



Let us consider the following Hamiltonian for a two-level system under the action of an external 
time-dependent field 

Htit) = Ha + Htit) = eas-fit)^ (1.1) 
and the corresponding Schrodinger equationP] 

= #i(£)*(£), (1.2) 

with ^ : R — > C 2 . Here /(£) is a function of time £ and e G R is a parameter representing half of 
the energy difference between the "free" (i.e., for / = 0) energy levels. The symbols a±, cr 2 and o" 3 
denote the Pauli matrices in their usual representations: 

ai = ( i I ) ' a2 = ( °i ~o ) and a3 = ( I -1 ) ' 

satisfying the commutation relations [<7i, = 2ia^, plus cyclic permutations. 

The "interaction Hamiltonian" Hi(t) := —f{t)a% represents a time-dependent external interac- 
tion coupled to the system inducing transitions between the two eigen-states of the free Hamiltonian 
Hq := eo"3. The situation where e is "small" characterizes the "large coupling domain" 0-0]. 

The system described above is certainly one of the simplest non-trivial time-depending quantum 
systems and the study of the solutions of ( |1.2p is of basic importance for many physical applications 
as, e.g., in quantum optics or in problems of quantum tunnelling. 

Equation ( |1.2| ) has been analysed by many authors in various approximations. In the wide 
literature on this subject we mention the pioneering work of Autler and Townes ||], where these 
authors studied the solutions of (|1.2|) for the case where, in our notation, /(£) = — 2/3 cos(co>£), 
f3 G R. Their work is exact but non-rigorous and involved a combination of the method of con- 
tinued fractions, for relating the coefficients the Fourier decomposition of the wave functions, with 
numerical analysis. No proof has been obtained that the continued fractions converge and further 
unjustified restrictions have been made in order to transform some transcendental equations into 
low order algebraic equations, which are then solved either exactly or, specially for strong fields, 
numerically. 

For a recent review on the mathematical theory of quantum systems submited to time- depending 
periodic and quasi-periodic perturbations see 0. For an introduction to the subject of "quantum 
chaos" and quantum stability, see |5| . See also || for a spectral analysis of the quasi-energy operator 
for two-level atoms in the quasi-periodic case. 

In we studied the system described by ( p..2| ) in the situation where / is a quasi-periodic 
function of time and a special sort of perturbative expansion (power series expansion in e) has been 
developed. Its main virtue is to be free of the so-called "secular terms", i.e., polynomials in £ that 
appear order by order in perturbation theory and that spoil the analysis of convergence of the series 
and the proofs of quasi-periodicity of the perturbative terms. Although we have not been able to 
prove convergence of our power series expansion in the general case where / is quasi-periodic it has 
been established that the coefficients of the expansion are indeed quasi-periodic functions of time. 

2 For simplicity we shall adopt here a system of units with h = 1. 
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One of the obstacles found in the attempt to prove convergence of the series is the presence 
of "small denominators". This typical feature of perturbative approximations for solutions of 
differential equations with quasi-periodic coefficients is well known as one of the main sources of 
problems in the mathematically precise treatment of such equations. 

On what concerns proofs of convergence it should, therefore, be expected that better results 
could be obtained if the function / were restricted to be periodic since, in this case, no problems 
with small denominators should afflict our expansions. 

However, the problem with small denominators is not the only problem to be faced in the 
perturbative expansion of [[[]]. In this paper we show how to circumvent the additional sources of 
difficulties and to finally establish convergence of our perturbative expansion for periodic /. 

By a time-independent unitary transformation, representing a rotation of ir/2 around the 2-axis, 
we may replace Hi(t) by 



H 2 (t) := (e— 2 / 4 ) H^t) (e^ 4 ) = er 1 + f(t)a 3 
and the Schrodinger equation becomes 

id t $(t) = H 2 (t)$(t), 

with 

$(£) : = e - i7R7a/4 #(t). 



;i.3) 



;i.4) 



;i.5) 



The theorem below, proven in presents the solution of the Schrodinger equation (|1.4|) in 
terms of particular solutions of a generalized Riccati equation. 



1.1 Theorem. Let f :R^R, f G C\R) and e G R and let g 
solution of the generalized Riccati equation 



Then, the function $ : 



G' - %G 2 - 2tfG + ie 2 
I 2 given by 



0. 



g G C 1 (1R), be a particular 

(1.6) 



U(t)$(0) = U(t, 0)$(0), 



where 



( R(t) (l + ig(0)S(t)) 



U(t) 



-ieR{t)S{t) 



-ieR(t) S(t) R(t) (l - i g(0) S(tj\ , 



with 



and 



R(t) 



cxp 



S(t) :- 



i / (f(r)+g(r))dT 



f R(t)~ 2 dr 
Jo 



is a solution of the Schrodinger equation (TA ) with initial value $(0) 



0+(O) 

Mo) 



gC 2 . 



(1.7) 

(1.8) 
(1.9) 

;i.io) 

□ 
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For a proof of Theorem \L.1[ see |l|. Let us briefly describe some of the ideas leading to Theorem 
|"TT| and to other results of PJ. As we saw in [I], the solutions of the Schrodinger equation ( |1.4j ) can 
be studied in terms of the solutions of a particular complex version of Hill's equation: 

tf(t) + (*7'(f) + e 2 + f(tf) <j>{t) = 0. (1.11) 

In fact, a simple computation (see [|l[]) shows that the components <f>± of $(i) satisfy precisely 

0'; + (+zf + e 2 + / 2 )0 + = 

(1.12) 

C + Hr + e 2 + / 2 )0_ = 



As a side remark we note that equations (|1 . 1 2|) are simpler and more convenient than the 
equations obtained by separating ip + and ip- from (|1.2|) : 



te 











;i.i3) 



This last pair of equations, mentioned (but not used) in f§], is mathematically less convenient 
because the coefficient /'/ / can be discontinuous and unbounded in typical cases as, for instance 
when f(t) = — 2(3 cos(cjt), the case analysed in ||. 

In [|l[ we attempted to solve (|1.11|) using the Ansatz 



</>(*) = exp l-i / (f(r)+g(r))dT 



;i.i4) 



It follows that g has to satisfy the generalized Riccati equation (|1.6|) and we tried to find solutions 
for g in terms of a power expansion in e like 



9(t) = qWj^fCn®, 



n=l 



where 



exp 



f(r)dr 



;i.l5) 



;i.l6) 



The heuristic idea behind the Ansatze (|l-14j) and ( 1.15 ) is the following. For e = a solution 
for ( [1.11] ) is given by exp ^— i J* /(r)drV Thus, in ( |1.14 ) and Ql.15 ) we are searching for solutions 
in terms of an "effective external field" of the form f + g, with g vanishing for e = 0. 



Notice that a solution of the form ( 1.14 ) leads to only one of the two independent solutions of 
the second order Hill's equation fll.ll ). The complete solution of the Schrodinger equation ( |1.4| ) in 
terms of solutions of the generalized Riccati equation (|1.6| ) is that described in Theorem |i~T . 



As mentioned above, perturbative solutions of quasi-periodically time-dependent systems are 
usually plagued by small denominators and by the presence of the so-called "secular terms". In 
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[p]] we discovered a particular way to eliminate completely the secular terms from the perturbative 
expansion of g (see Appendix ^) and we were able to show, under some special assumptions, that the 
coefficients c n (t) are all quasi-periodic functions. In [T|] we proved convergence of our perturbative 
solution in the somewhat trivial case where f(t) is a non-zero constant function. Unfortunately no 
conclusion could be drawn about the convergence of the perturbative expansion for g in the general 
case of quasi-periodic /. We conjectured, however, that our expansion is uniformly convergent at 
least in the situation where f(t) — M(f) is uniformly small. Here M(h) is the so-called "mean 
value" of an almost periodic function h, defined as (see, e.g. 0) 

M(h) := lim -J- / h(t)dt. (1.17) 



T^oo 2T 



T 



The technically central result of the present paper is the proof that, under suitable assumptions, 
the series ( P-.15Q converges absolutely and uniformly on 1 as a function of time for |e| small enough 
and / periodic. This is the content of Theorem B7TL Moreover, we show that the functions c n and, 
hence, g, have uniformly converging Fourier series representations. We use this fact together with 
the solution ( |1.8| ) to find the Floquet representation of the components <f>± of the wave function in 
terms of uniformly converging Fourier series representations. This is the content of Theorem |1.2| . 
Absolutely converging power series in e for the Fourier coefficients and for the secular frequency 
are also presented. 

We believe that the methods employed in this paper are also of importance for the general 
theory of Hill's equation. It would be of great interest to know whether the ideas described in |L| 
and here can be generalized and applied to a larger class of Hill's equations than those we studied 
so far. 



1.1 The Main Result 

On what concerns the solutions of the Schrodinger equation ( |1.4j ) the next theorem summarises our 
main results. 

1.2 Theorem. Let f be a real T^-periodic function of time (T w := 2-k/uo) whose Fourier decom- 
position 

f(t) = Y, F nC in "\ (1-18) 

with uo > 0, contains only a hnite number of terms, i.e., the set of integers {n G Z| F n ^ 0} is a 
finite set. Moreover, assume that F = 0. 

Consider the two following mutually exclusive conditions on f: 

I) M(q 2 ) ^ 0. 

II) M(q 2 ) = but M(Q X ) ^ 0, where 

Qi(t) := q{tf fq- 2 {r)dr. (1.19) 
Jo 

Then, for each f as above, satisfying condition I or II, there exists a constant K > (depending 
on the Fourier coefficients {F n , n6Z,n/0} and on u > 0) such that, for each e with \e\ < K, 
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there exist f2 e M and T^-periodic functions uf x and uf 2 such that the propagator U (t) of (T7i) can 
be written as 

U n (t) U 12 (t) \ ( U^t) U 12 (t) 
Fit) = | ] = [ J . (1.20) 

with 



Unit) U 22 {t) J \ -U 12 (t) U u (t) 



U u (t) = e- im u u (t) + e int u+(t), (1.21) 

U 12 (t) = e- im u u it) + e int uf 2 it). (1.22) 
The functions ufi and uf 2 have absolutely and uniformly converging Fourier expansions 

unit) = "£ U ti(n)e in "\ 



Moreover, under the same assumptions, and the Fourier coefficients Un( n ) and U^in) can be 
expressed in terms of absolutely converging power series on e. □ 



Remarks on Theorem 1.2 



1. Expressions ( 1.2 1|) and ( 1.22 ) represent the so-called "Floquet form" of the matrix elements 



Unit) and Uu(t). The frequency f2 is called the "secular frequency". 

In this paper we will assume that F = 0. Results on the almost resonant case F ^ 0, 
with Fq/lj satisfying some appropriated Diophantine conditions, will appear in a forthcoming 
publication Hl2| . 



3. The physically realistic condition that the Fourier decomposition of / contains only a finite 
number of terms can be weakened. The only condition we use is the fast decay for |m| — > oo 
of the Fourier coefficients Q m of the function q(t) (defined in (|1.16|) ), as found in Proposition 

4. The second equality in (|1.20|) is due to (|1.8|) . 

5. It is important to stress that conditions I and II are restrictions on the function / and not 
on the parameter e. 

6. Possibly there are other conditions beyond / and //which could be considered, but they have 
not been explored so far. They are relevant in some cases. Theorem |1.2j still does not provide 
a complete solution of ( |1.4j ) for all possible periodic functions /, but examples and some 
qualitative arguments show that the remaining cases are rather exceptional. For instance, 
for fit) = ifi cos(cjt) + ip 2 sin(c<jt) condition I covers all pairs ((pi, <p 2 ) G M 2 , except only the 
countable family of circles centered at the origin with radius x a u/2, a = 1, 2, . . ., where x a if 
the a-th zero of Jq in R + (J is the Bessel function of order zero). However, in these circles 
condition II is nowhere fulfilled. See the discussion in Section 0. 
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7. From the computational point of view the solution given by our method can be easily imple- 
mented in numerical programs and has been successfully tested, providing ways to study our 
two-level system for large times with controllable errors (due to the uniform convergence). 
Results on these numerical studies will be published elsewhere. 

8. Unitarity of U(t) for all t G R is a consequence of Dyson's expansion (see f.i. 0). 

9. Conditions I and II define, in principle, distinct solutions of the generalized Riccati equation 
(|1.6|) and, hence, of the Schrodinger equation ( |1.4| ). To fix a name we will call these solutions 
"classes" of solutions. 



1.2 Remarks on the Notation 

Let us make some remarks on the notation we use here and recall the notation used in Q. Given 
the Fourier representation?] 

/(*) = 2 (1-23) 

m6Z B 

of the quasi-periodic function /, we denote (as in [|lj) by u the vector of frequencies defined by 

( ujf G R B , if Fo = 

u := < ' . (1.24) 

[ (u f , F q )eR B+ \ if Fq^O, 

Since we assume that ujf G R+, the definition above says that all components of a> are always 
non-zero. Moreover, we denote 

B, if Fo = 
A := { ' . (1.25) 

B + 1, ifFo^O 

We denote vectors in 7h B (or R B ) by v and vectors in Z A (or R A ) by v. The symbol \n\ denotes the 
/ 1 (Z j4 ) norm of a vector n = (n 1; . . . , n^) G Z A : |n| := |rii| + • • - + \tia\- We use the symbol H for 
the identity matrix. Mat(n, C) is the set of all n x n matrices with complex entries. 

We denote by [^J the largest integer lower or equal to x G R. 

For m G Z we denote by <m> the following function: 

f \m\, for m ^ „„ N 
<m> := \ 1 " r . 1.26 

|^ 1, tor m = U 

For m G Z we denote by J m the Bessel function of first kind and order m. 

The symbol □ denotes end of statement and the symbol ■ denotes end of proof. 



3 For convenience we adopt here a different notation as that found in |IJ, where the Fourier decomposition of / 
was written as f(t) = ^ fme*W* 
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2 Some Previous Results 



In [|[ some results could be proven about the nature of some particular solutions of ( |1.6|) for the 
case where / is a quasi-periodic function subjected to some additional restrictions. These results 
are described in Theorem [2.1| . 

2.1 Theorem. Let f : R — ► R be quasi-periodic with 

fit) = 

and such that the sum above contains only a finite number of terms. Assume that the vector lj 



(defined in (jl.24) ) satisfies Diophantine conditions, i.e., assume the existence of constants A > 
and a > such that, for all n G 1, A , n^O, 

\n-u\ > A _1 |n| _CT . 

I. Assume that f satisfies the condition M(q 2 ) ^ 0. Then, there exists a formal power series 

oo 
n=l 

representing a particular solution of the generalized Riccati equation ( \l.dj) such that all coefficients 
c n can be chosen to be quasi-periodic and can be represented as 

c n (t) = ^e--', (2.2) 

where, for the Fourier coefficients , we have 

\C^\ < JC n e~^, 
where xo > is a constant and IC n > 0. 

II. Assume that f satisfies the conditions M(q 2 ) = and M(Qi) ^ 0, where Q\ is defined in ( |L lf\ ). 



Then, there exists a formal power series 

oo 

9(t) = q(t)J2 e n(t)e 2n , (2.3) 

n=l 

representing a particular solution of the generalized Riccati equation ( \1.6Q such that all coefficients 
e n can be chosen to be quasi-periodic and can be represented as 

e n (t) = Eg>e***, (2.4) 

me1 A 

where, for the Fourier coefficients , we have 

\E^\ < C p-x°H 

\ LJ m_ I — c n c ; 

where xo > is a constant and C n > 0. □ 
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There are other conditions beyond / and II which could be considered, but they have not been 
explored so far. See the discussion in Section [| 

The statements of this last theorem are not sufficient for proving convergence of the power series 
expansions in e for g. Unfortunately, as discussed in [jl]], the behavior for large n of the constants JC n 
and C n is apparently too bad to guarantee absolute convergence of the formal power series above. 

For the restricted case were / is periodic we will prove in the present paper stronger results 
(Theorem |3.1| below) than that implied by Theorem |2.1] . As we will see, these stronger results, in 
contrast, imply convergence of the e-power series for g (Theorem |3.3| below). 

Some of the more technical results of |L] have been obtained through the analysis of the Fourier 



coefficients of the functions c n and e n defined in Theorem ^TTj above. Specially important for us are 
the recursion relations found in [|TJ for the Fourier coefficients and Ej$ defined in (|2.2|) and 
(|2.4j) , respectively. Those recursion relations follow by imposing the generalized Riccati equation 
(|1.6| ) to the power expansions ( |2.1| ) and (2J3). In Appendix |A] we reproduce some of the main ideas 
of [l]] leading to a power series expansion for g free of secular terms and leading to the recursion 
relations below. 

It is important for our present purposes to reproduce those recursive relations here, what we 
shall do now. 



As in 0, let us denote by Qm the Fourier coefficients of the function q (defined in (|1.16|) ) 

q(t) = ^g m e im ^ (2.5) 



and by Q^m the Fourier coefficients of the function q 2 . For the Fourier coefficients of the functions 
c n we have found the following relations: 

Cfi> = a x Q m (2.6) 



/7(2) 



(~i{n) 



E 



(alQi > - QL 2 » ) 



E 



n ■ u 



1 



22.1: Xk2^^ \P — 




m— (n 1 +n 2 ) 



(2.7) 



Qn 



n-1 



for n > 3. 



(2.8) 



M(q 2 



Above m G Ij A , a\ = "\ K *„( . For the Fourier coefficients of the functions e n we have found the 



following relations. 



M(q 2 )' 



Qm+nQn ^ Q 



n ■ uj 



n (2) n (2) n (2) 

2iM(Q 1 ) ^ {n^ujir^-u) 



E 



(2.9) 
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Hi +212 #0 L 



Q, 



iM(Q 1 



(2) (2) 

Vn+tij +n 2 



n • a; 



n-1 



2iM{Qi 



+i-p) 



for n > 2. 



Above m G Z A , Qi is defined in ( |LT9|) and 

1 

K :-- 



2iM(Q 1 



n (2) n (2) n (2) 



—1 ' —2 



n-1 

p=i 

(n x + n 2 ) • a? 



(2.10) 



(2.11) 



The above expressions for the Fourier coefficients are somewhat complex but two important 
features can be distinguished. The first is the inevitable presence of "small denominators", rep- 
resented by the various factors of the form (n-u)~ l (with n ^ 0) appearing above. The second 
is the presence of convolution products (a consequence, lately, of the quadratic character of the 
generalized Riccati equation). The presence of the later is the additional source of complications 
mentioned before, for they also, together with the small denominators, contribute to spoil the decay 
of the Fourier coefficients needed to prove convergence of the e-expansions. 



3 The Recursive Relations in the Periodic Case 

In [U the recursion relations presented above have been used to prove inductively exponential 
bounds for the Fourier coefficients. As mentioned before two main difficulties have to be faced in 
this enterprise: the presence of "small denominators" and of convolution products in the recur- 
sion relations. Both are independently responsible for reducing the rate of decay of the Fourier 
coefficients at each induction step. 

Let us consider the origin of the "small denominators problem" in our recursion relations. 
It comes from the many factors of the form (n ■ (with n ^ 0) appearing in the recursion 

relations. In the case where / is a periodic function with frequency uj with F ^ 0, we have A = 2, 
n = (ni, n 2) G Z 2 and n • u — n^uj + n 2 F . On the other hand, in the case where / is a periodic 
function with frequency uj and with F = 0, we have A — 1, n = n<EZ and n • u = nuj. To avoid 
the quasi-resonant situation where nioj + U2Fq is small we will, as mentioned, consider in this paper 
the case where F Q = 0. 

For the Fourier coefficients of the functions c n , the recursive relations become 

Cg> = aiQ m , (3.1) 
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E 

n l+ n 27^° 



'n-1 



(ni + n 2 ) ■ uj \ ^— ' 



w m w n 



n-p) 



QmQ — ni 



Qm—(n\ 



+n 2 ) 



"«2 



<2, 



(2) 




0, 



n-1 



Above m G Z. 

For the Fourier coefficients of the functions e„ we have: 



for n > 3. 



(3.3) 



p(l) _ Qm-HnOrti Qm ST~^ 

"17^0 ni^O, n 2 #0 



n (2) n (2) n (2) 
Vni+n 2 V"l 

(nio;)(ri2a;) 



(3.4) 



E 

n x , n 2 6Z 



Qm—ni—n2 ~^ 



Q, 



iM{Qi 



n-1 



n 3 gZ 
n 3 ^0 



o (2) o (2) 

Vn 3 +ni+n 2 V"3 

n 3 u 



n-p) 
"2 



p=1 

(ni + n 2 )u; 



n-1 



2iM(Q 1 ) 



E E ^Mt 1 ^, 



for n > 2. 



(3.5) 



P=2 mez 



It is clear here that no "small divisors" appear in this case, since now > uj^ 1 for n^O. 

Hence, the convolution products are the only remaining factors eventually forcing the reduction of 
the decay rate of the Fourier coefficients at the successive induction steps. 

In the Section [| we will show how the effect of the convolution products can be taken under 
control. The result is expressed in the following three theorems. 



3.1 Theorem. Let f : 

with F = 0. 



be periodic with a finite Fourier decomposition as in ( FJa ) and 



Case I. Consider the Fourier coefficients Cffl satisfying the recursion relations (\3.1\ ), ( \3.2{ ) and 
( 3.5 ). Under the hypothesis that M(q 2 ) ^ we have 



-X\m\ 



(3.6) 



for all n G N, and all m G Z, where x > is a constant and Cm> is defined in (\1.2dj ). Above, the 
coefficients K n do not depend on m and satisfy the recursion relation 



C 2 



'n-1 \ /n-1 

E^n-p + [J2 K P K - 

+1-P 

,p=l / \p=2 



(3.7) 



with K\ = K2 = C\, where C\ and C2 are positive constants which can be chosen larger than or 
equal to 1. 
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Case II. Consider the Fourier coefficients Em'' satisfying the recursion relations ( |3.4|) and j \3.3j) . 
Under the hypothesis that M(q 2 ) = and M(Qi) ^ we have 



|£i n) | < KL 



-xM 



(3.1 



for all neN, and all m G Z, wiere x > is a constant. Above, the coefficients K' n do not depend 
on m and satisfy the recursion relation 



K' 



£-2 



n-1 



'n-1 



vP=2 



(3.9) 



with K[ = K' 2 = Si, where £\ and £2 are positive constants which can be chosen larger than or 
equal to 1. □ 



Theorem [3.1| will be proven in Section [|. The importance of the recursive definition of the 



constants K n given in ( |3.7| ) or ( |5T5| ) is expressed in the following theorem, which says that the 
constants K n grow at most exponentially with n. 



3.2 Theorem. Let the constants K n be defined through the recurrence relations ( |3.7p or p.Sj j. 
Then there exist constants K > and K > (depending eventually on f) such that K n < K K n 
for all neN. □ 



The proof of Theorem |3.2| is found in Appendix and makes interesting use of properties of 
the Catalan sequence. Theorems [3.1| and |3.2| have the following immediate corollary: 



3.3 Theorem. 

e £ C with lei 



The power series expansions in ( |2. 1\ ) and $2. 
K for all t G 



< 



and, hence, ( |2. 1\ ) and 



are absolutely convergent for all 
define particular solutions of the 



generalized Riccati equation (l^t) in cases I and II, respectively, of Theorem 3J. The function g 
can be expressed in terms of an absolutely and uniformly converging Fourier series whose coefficients 
can be expressed in terms of absolutely converging power series in e for all e e C with \e\ < K. □ 



Proof of Theorem |3.3| We prove the statement for case I. Case II is analogous. The first step is 
to determine the Fourier expansion of the function g, as given in ( |1.15|) , and to study some of their 
properties. One clearly has 

git) = ^G m e im ^, (3.10) 



with 



where 



G. 



Q(n) 



n=l 



^ ^ Qm-lCj 



(n) 



(3.11) 
(3.12) 
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Now and in future proofs we will make use of the following important lemma, whose proof is 
given in Appendix |C[ 



3.4 Lemma. For x > and m e Z define 

B(m) = B(m, x) ■= 

Then one has 

B{m) < B - 



-X(\m-n\+\n\) 



ne'i 



<Cm — n^> 2 <CrC^> 2 



~x\ m \ 



<m> 2 

for some constant Bo = Bq(x) > and for all m£Z. 



(3.13) 

(3.14) 
□ 



We have the following proposition: 
3.5 Proposition. For all x > there exists a constant C g = C g (x) > such that 



\G^\ < C g K, 



9 "<m> 2 

for all m & Z and all n e N. Consequently, for \e\ < K one has 

e -xH 

|G m | < C ff 



for some constant C (x, e) > and for all m G Z. 



(3.15) 



(3.16) 



□ 



Proof of Proposition 3.5 . Inserting (|3.6| ) and (4J.) into ( 3.12 ) we have, for any x > 

\G%>\ < K n Q B{m, x), (3.17) 
where B(m, x) is defined in (|3.13| ). Relation ( p,15|) follows now from Lemma |3~^ . ■ 



From this the rest of the proof of Theorem 3.3 follows immediately. 



The solutions for the generalized Riccati equation ( |1.6| ) mentioned in Theorem ^T3] are, through 



(|1.8|) , the main ingredient for the solution of the Schrodinger equation (|1.4|) . This will be further 
discussed in Section [5[ Now we have to prove Theorem |3~T . 



4 Inductive Bounds for the Fourier Coefficients 



In this section we will prove Theorem |3.1| in cases I and II. We will make use of the following propo- 
sition on the decay of the Fourier coefficients Q m and Qm of the functions q and q 2 , respectively. 
The proof of this proposition appears in Appendix |B[ 
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4.1 Proposition. Let f : R — > K be periodic and be represented by a finite Fourier series as in 
Then, for any constant x > there is a positive constant Q = Q(x) such that 

p-xH 

and 

P -X\m\ 

l«l < (4-2) 
for all m G Z, where <m> is defined in (\1.2t ). □ 



4.1 Case I 



In this section we will prove Theorem |3.1| in case I. Making use of Proposition [O] and of relations 
(|3.1|), (|3.2|) and (Of) we easily derive the following estimates: 



\C$\ < Q 



-xM 



l^ 2) l < 2uT 1 Q£ 



-x\m\ 



'n-l 



•x(H+|rn|) 



ni,n2eZ \p=l 



„— x|m-(ni+n2)| 



(4.3) 
(4.4) 



-X(l m l + i n l+™2|) 



n-l 



— V 771 

e v V irWi ir (n+1 ~ p) 
2I0 (2) I <™> 2 ^ ^ 1 1 

Z |V0 I ni£Z p=2 



(4.5) 



for n > 3. 



It follows from (P), from the definition of B{m) in (|3. 13|) and from Lemma |3.4| that 



|C7«| < Tu-iQ ( B(m) + ^^T, E 



-2x|ni 



< K 2 



-xH 



(4.6) 



|g ( 2)|<<m>> 2Z_, <<ni>> 

for some convenient choice of the constant K 2 . 

Now, we will use an induction argument to establish (|3.6|) for all n > 3. Let us assume that, 
for a given n G N, n > 3, one has 

g-xM 

|C£>| < K p — — - VmGZ, (4.7) 



for all p such that 1 < p < n — 1, for some convenient constants We will establish that this 
implies the same sort of bound for p = n. Notice, by taking K\ > Q, that relation ( |4.3[ ) guarantees 
(|4.7|) for p = 1 and that relation (O) guarantees the case p = 2. 
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From ( |4.5|) and from the induction hypothesis, 



'n-l 



-X(\ m -( n l+ n 2)\ + \ n l\ + \ n '2\) 



+ 



+ 



<m - (ni + n 2 )> 2 <Cni> 2 <n 2 » 2 



E 



' ni, n2£ 

p -x(|wi+n2|+|ni|+|n2|) 



|n( 2 ) I <m> 2 ^ , <ni + n 2 > 2 <Cni> 2 «n 2 » i 



Q e 



X\m\ 



'n-l 
vP=2 



E 



Now, 



E 



-x(l™l+™2| + |"l|+|«2|) 



ni, n2t 

are just finite constants and 



«ni + n 2 > 2 -Cni> 2 <n 2 >^ 



E 

rai,n2€Z 



,-x(l m -( n l+n2)| + |ni| + |n 2 |) 



e -2xlm| 



and 



xKI 



(4. 



E 



-2xl"i 



<m - (ni + n 2 )> 2 <ni> 2 <n 2 >^ 



-X(l(w-ni)-n2)|+|n 2 |) 



, <ni> 2 <(m - ni) - n 2 )> 2 <n 2 > 2 

niGZ 712 



E 



e -xlm| 

-£>(m - 

<ni> 2 



< 



B »E 



-x(l"-i|+l m -«i|) 



<ni> 2 <m - m»< 



< (Bo) 



-X\m\ 



(4.9) 



where we again used Lemma ^4 . 
Therefore, we conclude 



\C™\ < 



'n-l 



'n-l 



KpKn+l-p 



v p=2 



-X|m| 



(4.10) 



for two positive constants C a and C(,. Taking C 2 := max{C a , Cb, 1} relation fl3T7l) is proven with 
C 2 > 1. 

Notice that, without loss, we are allowed to choose K\ = K2 > 1 by choosing both equal to 
max{Ki, K2, 1}. ■ 
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4.2 Case II 



In this section we will prove Theorem 3J. in case II. From (|3.4j ) and (|3.5|) , from Proposition 4T and 
from the assumption ( |3.8|) we have 

O 2 ^ P -x(\m+m\+\ni\) 

\m < ^e 



«m + ni> 2 <Cni> 2 



+ 



Q4g-xH 



E 



-X(l n l+ n2 l + l n l 1 + 1^2 1) 



2 <m> 2 uj 2 \M(Qi)\ ^ , <m + n 2 » z <n x > z <n 2 > i 

' 1 ni, 7i2T 



(4.11) 



E 



Q 



-x(|m— m— nal+lnil+lnsl) 



+ 



ni,ri2GZ 

Q2g-x|m| 



<Cm — ni — n 2 3> 2 -Cni3> 2 <n 2 > 2 



-X(l™i+™2| + |ni| + |"2|) |7^| 



|M(Qi)| <m> 2 V<^i + ™2> 2 <rii> 2 <n 2 > 2 



Q e -X(l rl i+ n 2+ri3| + |"-i| + |ri2| + |ri3|) 

+ u <ni + n 2 + n 3 > 2 <ni> 2 «n 2 » 2 <n 3 > 2 



'n— 1 



n—p 



Qe 



X\™>\ 



E 



-2xM 



'n-1 



E W+i-p > forn>2. (4.12) 



vP=2 



Sums like 

E 



-X(|ni+n2| + |m| + |n 2 |) 



ni, ri2€Z 



<ni + n 2 > 2 <rii> 2 <n 2 >^ 



and 



-X(l n l+ rl 2+".3| + |".l| + |n.2| + |n3|) 



ni, ri2, ri3Sz 



<rii + n 2 + n 3 > 2 <rii> 2 <n 2 > 2 <n 3 >^ 



are just finite constants. By applying Lemma |3.4j we get 

e -xM 



\E&\ < €n 



(4.13) 



\Et ] \ < 



-X\m\ 



'n-l 



n— 1 



P=l / \p=2 



for n > 2, (4.14) 



where £ a , £b and £ c are constants. The rest of the proof follows the same steps of the proof of 
Theorem |3TI| in case I. ■ 



5 The Fourier Expansion for the Wave Function 



Now we return to the discussion of the solution ( |1.8| ) of the Schrodinger equation ( |1.4| ). Our 
intention is to find the Fourier expansion of the wave function $(£). 
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5.1 The Floquet Form of the Wave Function. The Fourier Decompo- 
sition and the Secular Frequency 

As explained in [[[]] and in Section [j], the components <fr± of the wave function $(£) are solutions of 
Hill's equation ( |1.12| ). For periodic / the classical theorem of Floquet (see e.g. [[UJ and |TlJ]) claims 
that there are particular solutions of equations like (|1.12p with the general form e lQt u(t), where 



u(t) is periodic with the same period of /. In order to preserve unitarity we must have Q e M. 
This form of the particular solutions is called the "Floquet form" and the frequencies Q are called 
"secular frequencies" . 

In this section we will recover the Floquet form of the wave function in terms of Fourier ex- 
pansions and we will find out expansions for the secular frequencies as converging power series 
expansions in e. 



According to the solution expressed in relation ( |1.7| ) and (|1.8|) , we have first to find out the 



Fourier expansion for the functions R and S defined in ( |1.9|) and (|1.10| ), respectively. 



We begin with the function R. The Fourier expansion of the function / + g is 

f(t) + g(t) = Q + J2( p n + G n {e)) e m "\ (5.1) 



nSZ 



where 
One has, 

with 



Q = Sl(e) := G (e). (5.2) 

(5.3) 



R(t) = e -^ (e) e~ m exp ( - H n e muJt ) 



Fra + Gn(£) , forn^O 
H„ = H„[r) : { , (5.4) 

0, for n = 



and 



7/(e) := iJ2 H m- (5-5) 

Notice that 7/(0) = 7/, where 7/ is defined in ( p.4j) . 

Since we are assuming that there are only finitely many non- vanishing coefficients F n , we have 



the following proposition as an obvious corollary of Proposition |3.5| : 



5.1 Proposition. For all x > and |e| small enough, there exists a constant Ch = Ch{Xi e) > 
such that 



e 



\H m \ < Ch—— 2 (5.6) 
for all m G Z. □ 
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Writing now the Fourier expansion of R(t) in the form 

-iUt \ r> „inu)t 



we find from 



-i7/ (e) 



1 + E E ^ ■ ■ ■ 



(5.7) 



for n = 0, 



-Rn = -Rn(e) — s 



1 ^") - ni,...,n p e 

p=l ni,...,n p € 



with 



jV p : = ^n a , 



(5.9) 



a=l 



for p > 1. 

In order to compute the Fourier expansion of S we have to compute first the Fourier expansion of 
R~ 2 . This is now an easy task, since the replacement R(t) — > R{t)~ 2 corresponds to the replacement 
(/ + 9) ~^ ~ 2(/ + g) and, hence, to if n — ► —2H n . We get 



J2(f) 



-2 _ 2ifii 



2) inujt 

& 1 



(5.10) 



with 



3 2i 7 /(e) 



3 2i7/(e) 



2 p+l 



1 + ^(p+l)! ^ 

p=l v 



• • • H n „H_ 



for n = 0, 



ni,..., n v 



2P^ 



2 ^ + E7-XT)! E H^-.-H^H^ , forn^O. 

(5.11) 

The following proposition will be used below. 

5.2 Proposition. For all \ > and |e| small enough, there exist constants Cr = Cr(Xj e) > and 
C R (-2) = C R (-2) (x, e) > such that 

e -xM 

\Rm\ < Cr—^ (5.12) 



l^- 2) l < c 



-X|m| 



R(- 2 ) 



for all m G Z. 



(5.13) 
□ 
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Proof of Proposition |5.2|. Using Proposition |5J] we have, for any p > 1 



m,...,n p €A 



< (C H ) P+1 Yl 6Xp ("^(Kl + • • • + Kl + \n-n 1 



ni,..., n p Gi 



Making repeated use of Lemma |3.4| on the right hand side of ( |5.14| ) we get 



^ H ni ■ ■ ■ H rip E, 



n-N z 



Inserting this into Q5.8|) gives (since _B > 1) 

I -Rn I < 



< 



B <n> 2 



3 | Im (7/(<0)|+Cff£o \ e -xN 



for all n G Z, as desired. The proof for Rn ^ is analogous. 



Assuming for a while 
we haveQ 



nu + 2Vt ^ for all nGZ, 



Arwjt 



with 



o n : — i 



p("2) 



and <t := — 5' n . 



f ' 
(5.14) 

(5.15) 



(5.16) 



(5.17) 
(5.18) 

(5.19) 



Assumption Q5.17 ) is actually a consequence of unitarity, as will be discussed in Section |5T2 . 
The following proposition is an elementary corollary of Proposition |5.2|: 



5.3 Proposition. For all \ > and \e\ small enough, there exists a constant C$ = C$(x, e) > 
such that 



\Sm.\ ^ C.t 



-x\m\ 



for all m e Z. 



(5.20) 

□ 



Writing 



17(f) 



C/n(t) M*) 
M*) ^22 (t) 



U u {t) U l2 {t) 



-U 12 (t) U n (t) 



For the case n = 0, (fj.lT]) says that f2 ^ 0. This must hold except for e = when = 0. 



(5.21) 
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we have for U\\ and Vyi- 



U n {t) 

u 12 (t) 



e~ lUt u u (t) + e lUt u+(t) 
e- int u- l2 {t) + e«* u+(t) 



with 



for 



with 



u u (t) 
r(t) 



= (l + ig(0)a o )r(t), u+(t) := ig(0)v(t), 
= —iea r(t), u i2(t) := —iev(t), 

^R n e inut and v{t) := ^ V n e i 



inwt 



ngZ 



(5.22) 
(5.23) 

(5.24) 

(5.25) 
(5.26) 



This provides the desired Floquet form for the components of the wave function $(£). We notice 
from the expressions above that the secular frequencies are ±f2. For Q we have the e-expansion 



(n) 
) 



n=l 



and for g(0), 

oo 

g (o) = £G ra = E en E G -- 

m&L n=l meZ 

Both converge absolutely for |e| < K, where K is mentioned in Theorem |3~2\ 



(5.27) 



(5.26 



As before, we have the following corollary of Propositions [572], |5.3| and Lemma |37J 



5.4 Proposition. For all x > and \e\ small enough, there exists a constant Cy = Cy(x, e) > 
such that 

e -xM 



iKnl < Cl 



for aii m € Z. 



(5.29) 
□ 



This last proposition closed the proof of Theorem 1.2. 



5.2 Remarks on the Unitarity of the Propagator 

The unitarity of the propagator U(t) means U(t)*U(t) = 11. After ( |5.21| ), this means 

\U n (t)\ 2 + \U 12 (t)\ 2 = 1. 



(5.30) 



Looking at relations (|5.22 ) and ( |5.23| ) two conclusions can be drawn from ( |5.30| ). The first is 
the following proposition: 
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5.5 Proposition. For eel and under the hypothesis leading to Qj5.22\ ) and ( \5.23{ ) one has OgR. 

□ 



The proof follows from the obvious observation that ( 5.30Q would be violated for \t\ large enough 
if Q had a non- vanishing imaginary part. Unfortunately a proof of this fact using directly the e- 
expansion of Q ( [5.271) is difficult and has not been found yet. 

The second conclusion is that ( 5.17| ) indeed holds. For, without this assumption there would 
be a term linear in t in (|5.18f) , violating (|5.30|) for large \t\. 

As in the case of Proposition |5.5| , no direct proof of this fact out of the e-expansion for Q (|5.27 ) 
has been found yet. The proof will probably follow the idea that is always smaller than 2u 
because Q is of order e and |e| has to be chosen small in order to provide convergence for the 
expansions. Analogously Q ^ because Q is analytic in e and, hence, has isolated zeros. If the 
analyticity domain must be small enough no zeros occur, except at e = 0. 



6 Discussion on the Classes of Solutions 



Let us now discuss some aspects of conditions I and II of Theorem |1.2j It is important to stress 
that these conditions are restrictions on the function / and not on the parameter e. 

As in (|B.1|) , let us write the Fourier decomposition of / as 

2.1 



f(t)=Y,f a e maU) \ (6.1) 



0=1 



with n a = —n 2 j- a +i and f a = f 2 j- a +i f° r an a with 1 < a < J. Comparing with Q1.18Q one has 
fa = F na , l<a<J. 

Hence, for F Q = and for fixed J and u, there are J independent complex coefficients f a and 
we can identify the parameter space M 2J with the set $j yU) of all possible functions / with a given 
J and uj. 

Condition M(q 2 ) = determines a (2 J — 1) or (2 J — 2)-dimensional subset of $j !tJ and there con- 
dition II applies. It is also on this subset that the more restrictive condition M(q 2 ) = M(Qi) = 
should hold, restricting the parameter space of / to a (2 J — 2), (2 J — 3) or (2J — 4)-dimensional sub- 
set. Hence, successive conditions like I and II would eventually exhaust completely the parameter 
space £ JiW . 

Conditions beyond I and II have not been yet analysed and many questions concerning the 
classes of solutions are still open. For instance, will further conditions like I and II really exhaust 
the parameter space of the functions /? Will the subtraction method of UJ and the convergence 
proofs of the present paper also work under these further conditions? What are the physically 
qualitative distinctions between the classes? Are these classes of solutions in some sense analytic 
continuations of each other? 

A distinction between class I and II may be pointed with the observation that in class I we have 
power expansions in e while in II we have power expansions in e 2 . Compare relations (|2.1|) and 



73T) of Theorem 2.1. 
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6.1 An Explicit Example 

To illustrate these ideas and point to some problems let us consider the important example where 
/ is given by 

fit) = <fi cos(ut) + if2 sin(a>i), (6.2) 

<p lt <p 2 e R. We have f(t) = fa~** + fa iwt with f x = (<Pi + m)/2, f 2 = J, J = 1, m = -1, 
ri2 = 1. Applying now ( p.5|) for this case with m = 0we get 

2p 



2^c 

CO 



(6.3) 



where cpo := w Vi + v 9 ! an d where Jo is the Bessel function of first kind and order zero. In this case 

7/ = vilu. 

Relation ( |6 .3| ) shows that condition I is not empty and that the locus in the ((fx, <£>2)-space of 
the condition M(q 2 ) = (necessary for condition II) is the countable family of circles centered at 
the origin with radius x a u/2, a = 1, 2, . . ., where x a if the a-th zero of Jo in R + . 

One shows analogously that 



Q, 



\\h\) m \ " 



and 



for all m G Z 

l(2) 



Q 



(2) 
m 



A 



(6.4) 



(6.5) 



For Q = the function Qi is periodic and we have in general 



M(Q, 



1,1 *-~< 



mgZ 



Q 



(2) 



m 



1 E 

r.i £-~< 



m=l 



( b (2) 



Q 



(2) 
m 



V 



rn 



Since |J m (x)| = |J_ m (x)| for all x G M, Vm G Z, it follows that \Qm\ 
for functions / like ( |6.2|) 

M(Q0 = 0. 



(2) 



(6.6) 



Vm G Z. Hence, 



(6.7) 



Therefore, condition II is nowhere fulfilled. For a complete solution of the problem for functions 
like ( |6.2| ), including the circles mentioned above, higher restrictions than that implied by condition 
II are necessary. 

6.2 A Second Example 

For functions / with J > 1 the situation leading to (|6.7| ) is not expected in general and condition 
II, and eventually others, may hold in non-empty regions of the parameter space of /. This can be 
seen in the following example with J = 2. Let us take 



ffl = h(t) + f2(t) 



(G.i 
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with 

fi(t) = he-^ + he^ (6.9) 

/ 2 (0 = f 2 e- i2wt + J 2 e i2 ^ (6.10) 
fi G C, i — 1, 2. We have g(t) = qi(t)q 2 (t), where 

gi (t) : = e* 7 * J^e inCl J n e ina; *, (6.11) 

q 2 (t) ■= e l ^J2e tnC2 J n (^\e m2 ^, (6.12) 



with 



l/il' 

It follows that 



e iCl = 7^, i = 1, 2. 



Q m = e^ + ^J2 el((m ' 2kKl+ ^ ( 6 - 13 ) 
Qg> = e 2i ^ + ^) ^e l «- 2fe )^ + ^)j m _ 2fe J fc . (6.14) 

From this we see (using J- n (x) = (—l) n J n (x)) that 

= ( _ l) m e -4 i(7/l+7/2 ) | e 2i( 7/l+7/2 ) ^ ( _ 1)feei (( m -2 fc )C 1+feC2 ) J m _ 2fc ^l^Tj J fc l . 

(6.15) 

The factor between brackets differs from Qm due to the presence of the factor (— l) k in the sum 

(2,) f2) 

over k e Z. Hence, we should rather expect |Qm | 7^ IQ-ml i n this case, what most likely implies 
M(Qi) 7^ for M(q 2 ) = 0, leading to a non-empty condition II. 
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Appendices 



A Short Description of the Strategy Followed in [jtj 



For convenience of the reader we reproduce the main steps of the strategy developed in U for 
finding a power series solution of the generalized Riccati equation (|1.6|) without secular terms. 

As discussed in Section [I], a natural proposal is to express g, a particular solution of (|1.6|), as a 
formal power expansion on e which vanishes at e = 0. For convenience, we write this expansion as 
in ( |1.15| ) where q(t) is defined in ( |1.16| ). This would give the desired solution, provided the infinite 
sum converges. Inserting ( 1.15 ) into ( |1.6| ) leads to 



n-l 



p=l 



n=l 



Assuming that the coefficients vanish order by order we conclude 

(qci)' - 2ifqc 1 = 0, 
(qc 2 )' - iq 2 cl - 2ifqc 2 + i = 0, 

n-l 

(qc n )' - i ^2 (f c P c n-p ~ 2ifqc n = 0, n > 3. 
P =i 



The solutions of ( |A.2| )-( pO| ) are 

ociq(t), 



Cl (t) 



Cn{t) = q(t) 



i f (alq(t') 2 -q(t')- 2 )dt' + a 2 
Jo 

V / cp(0cn-p(0 dt')+a 

p=l J ° J 



for n > 3, 



(A.l) 



(A.2) 
(A3) 

(A.4) 



(A.5) 
(A.6) 

(A.7) 



where the a n 's above, n — 1, 2, ... , are arbitrary integration constants. 

The key idea is to fix the integration constants in such a way as to eliminate the constant 
terms from the integrands in (|A.6| ) and ( |A.7|) . The remaining terms involve sums of exponentials 
like e inujt , n^O, which do not develop secular terms when integrated, in contrast to the constant 
terms. For instance, fixing at\ such that M(afq 2 — q~ 2 ) = 0, that means, a\ = M \q~ 2 ) / M \q 2 ) , 
prevents secular terms in ( |A.6 ). 
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As shown in this procedure can be implemented in all orders, fixing all constants oti and 
preventing secular terms in all functions c n (t). In case I, relations fl2.6p -(2T%) represent precisely 
relations ( |A.5|) -( |A77 ) in Fourier space with the integration constants fixed as explained above. Case 
II is analogous. 



B The Decay of the Fourier Coefficients of q and q 2 

To prove our main results on the Fourier coefficients of the functions c n and e n we have to establish 
some results on the decay of the Fourier coefficients of q and q 2 . 

We write the Fourier series ( |1.18| ) of / in the form[] 



fit) = ^P n e 



inujt 

> 



with F n = F_ n , since / is real. In order to simplify our analysis we will consider here the case 
where the sum above is a finite sum. This situation is physically more realistic anyway. 

By assumption, the set of integers {n G Z| F n ^ 0} is a finite set and, by the condition that / 
is real and Fq = 0, it contains an even number of elements, say 2 J with J > 1. Let us write this 
set of integers as {ni, . . . , n 2 j} and write 

2J 

/(t) = X)/« einohrt ' ( R1 ) 

a=l 

with the convention that n a = — n 2 j_ a+ i, for all 1 < a < J, with f a = F na . Clearly /„ = faj- a +\, 
1 < a < J. 

A simple computation (see 0) now shows that q has a Fourier decomposition as in (|2.5|) with 

oo 2J 



pi,...,p 2 j=0 a=l 



fa 



Pa 



(B.2) 



where 



and where 



2.1 



P = P{p 1 ,...,p 2J ,n 1 ,...,n 2J ) := ^p t 6Z, (B.3) 



6=1 



2J f 

v := i^J!-. (B.4) 



a=l 



As one easily sees, 7/ € R. Above 5 (P, m) is the Kronecker delta: 

, \ / 1, if P = m, 
5(P ' m) := I 0, else. 



5 As above, here we adopt Fq = 0. 
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Since the function q 2 is obtained from q by replacing / — > 2f we have from 

(B.5) 



11 r ' '2/ e N "' 



Pa! \n a uJ 



Pl,—,P2J=0 a=l LJ 

where Q$ are the Fourier coefficients of q 2 . The coefficients Q m and Qm"* can also be expressed in 
terms of Bessel functions of the first kind and integer order. See Section |6] for some examples. 

As in |TJ, define 

fa 



09 := max 

l<a<2J 

and 

2J 



n n u) 



n b \. 



6=1 



Notice that, since the n&'s are fixed by the choice of /, Af is non-zero. 
The following important bounds have been proven in Appendix D: 

|Q m | < (2j e ( 2J -^) — (1--— ^- ) , (B.6) 

lv 1 _ v ; [AT" 1 ml! V TAT- 1 ml +1/ v ; 



and 



[A^Hl 1 V [A^H 

for all m with ["A/" _1 |m|] + 1 > 2(p. Above \x\ is the lowest integer larger than or equal to x. 

In we derived from ( |B.6[ ) a simple exponential bound for \Q m \, namely, 

\Q m \ < Qe- X|ro| , (B.8) 

where Q and x are some positive constants. For the purposes of this paper a sharper bound than 
is needed and we have to study relation ( B.6| ) more carefully. The result is expressed in 



Proposition |4.1| whose proof we present now. 

Proof of Proposition |4.1| . Let us consider first the coefficients Q m . Due to the dominating factor 
|A/" -1 |m|] !, one has 

,. «m» 2 ^-vn 

hm j — ; — tt—i — — — — = 0. 



| m |^oo e -xM [tV-i 



for any constant x > 0. Hence, one can choose a constant M\ > depending on x such that 

mPA^Ml P -x\m\ 



\J\f- 1 \m\]\ <m> 2 

for all m G Z. Therefore, there exists a positive constant Qi > (depending on x) such that 
\Qm\<Qi <m>- 2 e"*H for all m G Z. For we proceed in the same way and get the bound 
\Qm\< 0,2 <m>~ 2 e~*H for all m G Z. In Q and (p|) we adopt Q = max{Qi, Q 2 }. ■ 



Remark. The proof of Proposition [4.1| shows that we have also sharper bounds like 

e -xH 

\Qm\ < Qk—-—j: 

for any k G N. For the purposes of the present paper it was enough to consider k = 2. 
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C Bounds on Convolutions 



Here we will prove Lemma [3.4| . Consider for x > an d ui G Z 

p -x(|w-n|+|n|) 

B(m) = B(m, x) : = £ rj^Ta- (CI) 

First notice that B(m) = B(—m) for all m G Z. Choosing _B to be such that 

e -2xM 



7 ««>" 

the statement of the lemma becomes trivially correct for m — 0. Hence, it is enough to consider 
the case where m > 0. 

In flC.l| ), the sum over all n G N can be split into three sums: 

-1 2 X n ™ i °° -2 X n 

B(m) = e-* m y - +e -* m Y - + e* m V — - (C.2) 

^— ' (to - n) 2 n 2 ^— ' <m - n> 2 <n> 2 ^— ' (to - n)% 2 

n=— oo v ' n=0 n=m+l v ' 

In the first sum above we perform the change of variables n — > —n and in the third sum we perform 
the change of variables n — » n + m. The result is 

/ 00 -2 X n m 1 \ 

\ ^ (m + n) 2 n 2 ^ <m - n> 2 <n> 2 / 

\ n=l v ' n=0 / 

Now we will study separately each of the sums in (|C.3| ). Since for n > 1 one has m + n> <Cto^> 
one has for the first sum 

V — < 1 (C.4) 

^— ' (m + n) 2 n 2 <^jri^> 2 

n=l v ' 



where B\ := 



r? 2 



n=l 

The second sum in ( |C.3| ) is a little more involving. We have 

m j L m /2J ^ m ^ 

E <m - n> 2 <n> 2 ~~ E <m - n> 2 <n> 2 + E <m - n> 2 <n> 2 

n=0 n=0 n=[m/2j+l 

For the first sum in the right hand side of ( |C.5| ) we have -Cm — ri^> > in — n > m— [m/2\ > to/2. 
For the second sum in the right hand side of (|C.5|) we have n > [m/2\ +1 > m/2. Hence, for 
m > 0, 



™ 1 < /2_^ 2 

n=0 x 



Lm/2J 



E^^+ E 



<n> 2 ^ <gjn — n;» 2 

n=0 n=|m/2j + l 



< 2 f— Ve^ ( c - 6 ) 

x 7 n=0 



2S 



Therefore, choosing 



the lemma is proven. 



X 

B = 2B X + 8 V 



n=0 



(C.7) 



The proof of this lemma has the following proposition as corollary, generalizing Lemma [374 
C.l Proposition. For x > 0, k e N, k > 2, let 

Then, there exists a constant B 0t depending eventually on k, such that 



B k {m) < B ,, 



<m> A 



for all m G Z. 



(C.9) 

□ 



D Catalan Numbers. Bounds on the Constants K 



n 



Here we will prove Theorem |37J|. Let us start recalling that we have chosen K x = K 2 = C x for some 
constant C x which, in turn, can be chosen without loss to be larger than or equal to 1. The proof 
of Theorem |3.2| will be presented on four steps. 

Step 1. In this step we show that the sequence K n , defined in ( p.7|) , is an increasing sequence. 

First notice that K 3 = C 2 (2K X K 2 + (K 2 ) 2 ) = 3C 2 (K 2 ) 2 . Since K x = K 2 > 1 and C 2 > 1, we 
have Ki = K 2 < K 3 . 

Let us now suppose that 

K x = K 2 <K 3 <---<K n (D.l) 
for some n > 3. We will show that K n+X > K n . We have 

i n n—X Ti—l 

^ KpK n _ p+ i + ^ ^ K p K n - p+2 — ^ K p K n _p — K p K n _ p+ i 



K n +i — K n — C 2 



p=2 



p=l 



p=2 



71-1 



2K x K n + 2^ KpK n ^ p+2 — 2^ KpK n ^p 

p=2 p=l 



= C 2 [2K x K n + (if 2 K n - K n _ 2 Ki) + (jr 3 - + • • • + (K n - K n _ 2 )K 2 ] 

= C 2 [2K x K n + (K n - K n „ 2 )K x + (tf 3 - KO^.! + ...+ (K n - K n ^ 2 )K 2 ] , 

where in the last equality we used K x = K 2 . Now, from hypothesis ( |D.1| ) we conclude that 
K n+X > K n , thus proving that K n is an increasing sequence. 
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Step 2. Here we show that the sequence K n defined in ( |3.7|) satisfies 

n-l 

K n < 3C 2 J2K p K n . p+1 (D.2) 

p=2 

for all n > 3. 

We have already shown that K 3 = 3C 2 (K 2 ) 2 . Hence, ( |L).2| ) is obeyed for n = 3. 



Assume now that ( p.2|) is satisfied for all K p with p & {1, . . . ,n — 1}, for some n > 4. We will 
show that it is also satisfied for i^ n . In fact, we have from ( |3.7| ) 

K n = C 2 [K x K n . x + K 2 {K n _ 2 + K n _ x ) + K 3 (K n „ 3 + K n _ 2 ) + • • • + + K 2 )] . (D.3) 

From this and from the fact proven in step 1 that the sequence K n is increasing, it follows that 

K n < C 2 [K x K n _ x + 2 (K 2 K n . x + K 3 K n „ 2 + ■■■ + K n _ x K 2 )} (D.4) 

Now, using the obvious relation 

K x K n _ x = K 2 K n . x < (K 2 K n . x + K 3 K n _ 2 + ■■■ + K n _ x K 2 ) 

we get finally from ( p.4|) 

n-l 

K n < 3C 2 [K 2 K n _ x + K 3 K n _ 2 + --- + K n _ x K 2 ] = 3C 2 J2 K P K n- P +i, (D.5) 

p=2 



thus proving 

Step 3. Here we will prove the following statement. Let L n be defined as the sequence such that 
L x = L 2 = K x = K 2 = C x and 

n-l 

L n = 3C 2 ^ ^ LpL n - p + x . (D.6) 

p=2 

Then, one has 

K n < L n , Vn G N. (D.7) 

First notice that K 3 = 3C 2 {K x f = 3C 2 {L X ) 2 = L 3 . Hence, (|IX7| ) is valid for n G {1, 2, 3}. Now 
suppose K p < L p for all p G {1, . . . , n — 1} for some n > 4. One has from (p.2|) 

n— 1 n— 1 

K n < 3C 2 KpK n ^ p+x < 3C 2 L p L n - p+x = L n , (D.8) 

p=2 p=2 



thus proving ( p.7|) . 

Step 4. Consider the sequence c n defined as follows: c x = c 2 = 1 and 



n-l 

5^c p c n _ p+ i (D.9) 

p=2 
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for n > 3. The so defined numbers c n are called "Catalan numbers", after the mathematician 
Eugene C. Catalan. The Catalan numbers arise in several combinatorial problems (for a historical 
account with proofs, see 0) and can be expressed in a closed form as 

- - (2 "" 4)! n>2. (D.10) 



h - l)\(n - 2) 



(see, f.i, or ||). Using Stirling's formula we get the following asymptotic behaviour for the 
Catalan numbers: 

1 4 n 

C " W 16^ ^' " larga (D11) 

The existence of a connection between the Catalan numbers and the sequence L n defined above 
is evident. Two distinctions are the factor 3C 2 appearing in ( p.6| ) and the fact that Li = L 2 = C\ 
is not necessarily equal to 1. Nevertheless, using the definition of the Catalan numbers in ( p.9|) , it 
is easy to prove the following closed expression for the numbers L n : 



L n = (d)"- 1 (3C 2 r 2 (2 " 4) * , , n>2. (D.12) 

(n — l)l(n — 2J! 

We omit the proof here. Hence, the following asymptotic behaviour can be established: 

From the inequality K n < L n , proven in step 3, it follows that K n < Kq(12CiC2) u for some 
constant Kq > 0, for all nGN. Theorem |3.2j is now proven. ■ 
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